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[T I I Abstract 

In this note, we propose an expression for the eigenvectors and scalar products 
for a class of spin chains with long-range interaction and su{2) symmetry. This 
class includes the Inozemtsev spin chain as well as the BDS spin chain, which 
is a reduction of the one-dimensional Hubbard model at half-filling to the spin 
f^ ; sector. The proposal is valid for large spin chains and is based on the construction 

QJ • of the monodromy matrix using the Dunkl operators. For the Inozemtsev model 



Introduction: Integrability is a powerful tool in theoretical physics. The 



these operators are known explicitly. This construction gives in particular the 
eigenvectors of (an operator closely related to) the dilatation operator of the AA = 4 
gauge theory in the su{2) sector up to three-loop order, as well as their scalar 
(Nj . products. We suggest how this will affect the expression for the quasi classical 

'nI" I limit of the three-point functions obtained by I. Rostov and how to include the 

r^ ■ all-loop interaction. 

cn 

o 

p^ . most spectacular recent achievement it underlies is to determine the spectrum of 

the A/" = 4 SYM theory for an arbitrary coupling constant and to make very precise 
quantitative predictions confirming the validity the AdS/CFT correspondence [1]. 

'k^ • For recent reviews on this subject see [U [3]. The correlation functions [H El El 

^ , El E] ; the amplitudes [9] and the Wilson loops [TOl E] seem also to be within the 

reach of integrability methods. At weak coupling, the problem of computing the 
three-point function of single-trace operators at tree level in the su{2) sector was 
reformulated [12l |H E] in the language of spin chains in terms of scalar products the 
XXX spin chain. The scalar products are expressible as determinants [131 [HI [15] . 
The determinantal formula initially proposed for single-chain scalar products was 
generalized |16j to the overlaps involving several spin chains. Several special cases 
were considered [5l [6] involving at least one operator which is close to BPS. The 
main difficulty in taking the classical limit, which is necessary to compare with the 
string predictions, resides in taking continuum limit of large Slavnov determinants. 
This difficulty was overcome very recently by I. Kostov [T7], who succeeded to 
obtain a factorization formula for the Slavnov determinants. He obtained the 
continuum limit of the correlator of three operators with three sets of su(2) charges. 



At higher loop order, the dilatation operator is mapped to a spin chain with long 
range interaction [TH] . In [TO] , it was matched with the Inozemtsev model [2D] up to 
three loop order. Another related chain which reproduces the dilatation operator 
up to three loops is the BDS spin chain [21j which can be be obtained by reduction 
[24) of the one-dimensional Hubbard model at half- filling to the spin sector. Both 
these models belong to a class of deformations of the XXX model considered in [25] . 
These models were considered out of reach of the algebraic Bethe ansatz techniques, 
which allows to build the eigenvectors of the of the Hamiltonian, because the 
monodromy matrix was not known. The sole exception was the Inozemtsev model 
at infinite length [19] . whose monodromy matrix can be obtained by analogy with 
the one for the Haldane-Shastry model [26 [I271[25] . The construction is based on an 
(iso)morphism with the inhomogeneous XXX model. Here, we work out the details 
for the Inozemtsev spin chain and we postulate that an identical construction 
underlies all the models which are constructed from the boost charges in [25] . 
including the BDS model. The Hamiltonian obtained from the monodromy matrix 
at finite length has a periodicity defect but, on the basis of the perturbative results 
obtained by Gromov and Vieira for the periodic finite chains [8j, we think that it 
is possible to modify the construction to include perfectly periodic chains. It is 
not clear whether the lack of periodicity of the Hamiltonian really is a drawback 
in the context of the computations of the correlation functions in the A/" = 4 SYM 
theory, since the procedure of splitting and joining of the chains effectively breaks 
the periodicity. We certainly need to understand the origin of the cancellation 
between the contact terms necessary to restore the periodicity and the Hamiltonian 
insertions at the splitting/joining points which was observed at one loop in [8j. This 
work was triggered by discussions with N. Gromov and P. Vieira and by their recent 
results ^ on the eigenvectors of the dilatation operator and the scalar products 
at two and three loops. 

The Inozemtsev spin chain: The Inozemtsev spin chain [20] is an integrable 
long-range interacting spin chain with Hamiltonian 

L L-l 



H = Y.Y. T'L,.l.in){\ - P,,,+n) , (1) 



j=l n=l 

where Vi^^/f^{z) is the Weierstrass function with periods L and iir/n and Pj,j+n 
permutes the spins at sites j, j + n. In [19j it has been shown that around the 
limit K — >• oo this Hamiltonian matches the dilatation operator of the AA = 4 gauge 
theory in the su{2) sector up to three-loop order. In fact, the dilatation operator 
is a combination of the Hamiltonian ([1]) and a higher conserved quantity, but this 
aspect is irrelevant to the construction of the eigenvectors. The full solution of the 
Hamiltonian ([1]) is rather involved and its integrable structure was not yet fully 
explored. The limit k = corresponds to the Haldane-Shastry spin chain, whose 
monodromy matrix was constructed in [27J. A systematic way to construct the 
integrals of motion was given in [^. Another limit which is tractable is L — )• oo, 
when the interaction strength becomes proportional to 1/sinh nn. In this case, 
the conserved quantities can be obtained [2^ from those of the Haldane-Shastry 
model by the exchange of the real and imaginary periods L — )■ m/ k,. The Yangian 
generators are also obtainable by this procedure. The two limiting cases mentioned 
above, Haldane-Shastry and infinite length Inozemtsev model share the particular 



feature that the Yangian is a symmetry of the Hamiltonian. This can be understood 
easily from the fact that the conserved quantities are determined [25] by (a special 
limit of) the quantum determinant, which commutes with the monodromy matrix. 
When K — )• cx), we retrieve the conserved quantities of the XXX spin chain. Here, 
we are interested in exploring the first few corrections to the XXX limit. The 
expansion parameter will be the gauge theory coupling constant g'^, related to k 
by the relation [19] 

e-2« =g^- 3/ + 0(/) . (2) 

The strategy of the construction is as follow: we use the mondromy matrix 
of the Haldane-Shastry spin chain and then we continue it analytically in the 
positions of the spins, Zj = e'^'"^^'^ — >• e'^^'^. The algebra between the generators is 
not affected by this change (but the spectrum of the operators will be) . The key 
elements which allow to incorporate the long-range nature of the interaction |27j 
are a set of mutually commuting operators known as the Dunkl operators {22} [23] 

[di,d,]=0. (3) 

In the case of Inozemtsev spin chain, they have the following expression 

j;j>i j;j<i j;J9^i j\i<i 

where 0jj = Zi/{zi — Zj) and Zj = e^^'^. The operators Kij permute the coordinates 
KijZj = ZiKij. Under the coordinate permutations, the Dunkl operator behave as 

[Ki,i+i,dk]=0 if k^i,i + l, (5) 

Ki^i+idi — dj+iKj,j+i = 1 . 

The Dunkl operators of the Haldane-Shastry and Inozemtsev spin chains are par- 
ticular in the sense that they are linear in the permutations Kij . The fact that they 
commute is based on the following identity satisfied by the functions Qij defined 
above 

Qij&ik = Qij&jk + Qik&kj (6) 

for any three indices i,j, k. If we relax this condition by taking Qij = 1/(1— e(j— i)), 
with e{n) some function, the commutation relations ([3]) will not satisfied and we 
will have to add to dj's higher odd powers in the permutations to compensate for 
this defect 

di = di+ Y. Vi;3kiKijKikKu + 0{K'>). (7) 

The coefficients rj i-ji^i will be determined in terms of Qij by the condition of van- 
ishing of the terms quadratic in Kij in the commutators. We believe that this 
procedure will allow to determine recursively all the coefficients in the expansion 
([7]) and to explicitly build the integrable long-range spin chains with multi-spin 
interaction we describe below. 

Given the Dunkl operators, one can easily construct the monodromy matrix 
|27) . The construction in the rest of this section is rather general, and at this point 
we do not have to specify the precise details of the Dunkl operators, except that 



they are built from the coordinates permutations Kij, as explained above. We 
define then 

uu, . .(ft („)) , w = n (i + „Jg-v2 ) • («) 

where a stands for the auxiliary space. The projection operator tv acts in the 
following way: the coordinate permutations Kij are brought to the right of the ex- 
pression and then they are replaced with the spin permutations Pij . This amounts 
to working on (wave) functions which are symmetric by simultaneous permutations 
of spins and coordinates, and then freezing the coordinates to some specific values. 
It was shown in [27j that T(n) obeys the Yang-Baxter equation with the rational 
R matrix R{u) = u + iP. The proof goes in two steps: first, one checks that 
the matrix ([SD without the projection obeys the Yang Baxter equation, which is 
obvious. Second, it can be shown using the commutation relations ([5]), see for ex- 
ample eq. (2.24) in reference [27], that the unprojected monodromy matrix T(u) 
preserves the space of function symmetric by simultaneous permutations of spins 
and coordinates, so that 

7v{f a{u)f a' (v)) = 7r(f„(n))7r(f„, (t;)) . (9) 

In other words, the BGHP projection [27j is a morphism. In virtue of this relation, 
the projected matrix obeys the Yang Baxter equation as well. For the Inozemtsev 
spin chain, the matrix ([8]) is defined when L — t- oo, and is not known how to 
build the monodromy matrix for a periodic finite chain. However, since we are 
interested in expanding the Hamiltonian in powers of g'^, according to ([2]), on a 
finite chain we can truncate the sum ^ to sites from 1 to L. The price to pay is 
that the resulting Hamiltonian, while remaining integrable, will not be perfectly 
periodic anymore and some terms connecting the first few and last few sites will 
be missing. This Hamiltonian will be slightly different from the one considered in 
[S], which probably explains the difference with the contact terms which appear 
in that reference. Here we will concentrate on (long) spin chains with the integral 
structure coming from ([S]), and we hope to elaborate more on the finite chains and 
periodic boundary conditions elsewhere [31]. 

The normalization of the matrix ^ is chosen such that at dj = we retrieve 
the standard expression of the monodromy matrix of the XXX modeU . Prom the 
definition of the projection tt it follows that the symmetric combinations of Dunkl 
operators, for example the symmetric sums, have the property 

7r(F(d) j;<) =7r(F(d))7r(j;dr) , (10) 

i i 

where F(d) is an arbitrary function of dj's. We can therefore define as mutually 
commuting quantities the symmetric sums 



Pn = '^[^d") , [p„,Pm]=0. (11) 



^ The i in front of dj insures the right ratio between Pqj and dj , cf. [27] . 



They cannot be used to generate the Hamiltonians of the model, because they have 
the the same value on all the states |27l [2H] . This means that their eigenvalues can 
be determined on any state, for example on |r2) = | ff . . . f). 

Prom the expression ([8|) its is straightforward to obtain the Yangian generators 
by considering the coefficients of the expansion around n = oo, 



Qt 



ah 



Z,Ef 



(12) 



Qt = EJ>^ErEf + EJEf7^id 



'J J ■> 



(13) 



Tpab Tpba 



where E'^ are the elementary generators of gl{2) so that Pij 

Let us mention here that the construction of the monodromy matrix based on 
Dunkl operators can be done completely similarly for the gl{n) case. The boundary 
case works as well, as explained in |30) . There, two type of boundary reflection 
matrices were considered, one preserving the su{2) symmetry, the other which 
breaks it. 

The monodromy matrix ([8]) strongly resembles the monodromy matrix of the 
inhomogeneous XXX model defined as 



To(u;i 



n 



1 + 



^-^ la 



U ■ 



i/2 



(14) 



It would be tempting to identify 9j with the eigenvalue of dj, but this cannot be 
done because of the projection tt. However, there is a way to relate the two models. 
Let us first use the monodromy matrix acting on the reference state |0). Since 
the Dunkl operators mutually commute and commute with the spin operators, we 
can formally expand the matrix T{u) in powers of dj's. Up to trivial powers of i, 
the coefficients of the expansion are the same as for the expansion of Tq{u;6) in 
powers of 6'j's, 



num) = Y: E E 






n ji<...<jn ki,...,k„ 



kA...kJ 



a;/...9;;ro(n;0) 



TV 



e=o 



or 



T{u)\n)=VeToiu;e)\n) 
where the theta operator defined as 



Vf, 



E E 

" jl<---<jr, 



E 



ofcl- 



k ' 



Si,. 



e=o 






l^n 



Jl jn 



(15) 



(16) 



•■ji^n 



c^:::ri^)--K^---C^ 1^)' 



(17) 



bears strong resemblances to the theta derivative defined in [8j at two-loop and 
three-loop ordecl- From the property ([10]) of the projection we deduce that, if G{6) 



^In a previous version of this article, we have stated that the two differential operators are the same. 
In fact, we became aware that there are some subtle differences between the present definition and the 
one in [5], concerning in particular the terms which are odd in the derivatives, see below. We thank N. 
Gromov for explaining their definition of the differential operator. 



is a symmetric function of 9, we obtain |f| 

Vg[Fie) G{e)]\e=o = Ve[F{e)]\g=o Ve,[Gie')]\g,=o ■ 



in 



The action of the theta operator on any symmetric function of the impurities can 
be obtained by working out to the action of Dg on the basis of symmetric functions 
Pni • • • Prik with pn = J2j ^j" ^nd po = L. 

Let us define a more general operator ^g, via 



Tin) 



% Toiu; 



\e=o 



(19) 



E E E 

n ji<...<jnkl,...,kn 



ofcl- 



hi. 



k I 



5, 



fcl 



Jl 



..dl^Toiu; 



e=o 



TV d 
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Jn 



so that we can handle the monodromy matrix without acting on a particular state 
on the right. The coefficients of the differential operator are now spin operators, 
and they cannot be replaced by numbers anymore. Due to the morphism property 
([9]) of the projection tt, we deduce that the operator &g is also a morphism, which 
we call the Bernard-Gaudin-Haldane-Pasquier (BHGP) morphism 



Taiu)Ta'iv) = % n,a{u- 6) % T^.a'{u- 



6»=0 



= %\T^^aiu-,e)T^^a'iv-M (20) 



Let us emphasize that the purely differential operator T)g does not possess this 
property, rather 



Pe[ro,a(^;^)roy(^;^)] 



= Vg Toa{u; 9) Vg Tq „' [u] 6) + cross terms . (21) 
e=o e=o ' 61=0 



In (1201) the cross terms are taken into account by the commutator of tt (d-^ ■ ■ -d- 

with TQa'{u;6). Upon acting on the vacuum state we retrieve the original theta 
operator action 



%[...]{u-e) \n) = Vg[...]iu;0) \n) 



e=o 



(22) 



where the dots stand for any product of elements of the monodromy matrix. 

T^vo and three loops from the Inozemtsev model: Let us give an example 
how to perturbatively compute the operators df and dj dj and their projection on 
the vacuum state |0) using their expression in the Inozemtsev spin chain (j3|) . The 
connection with the dilatation operator of the AA = 4 theory works up to three 
loops [12], so we are going to expands to this order. From their definition, it is 
clear that the coefficients Qij are invariant by translation, 0jj = Qi^nj+n- The 
projection tt amounts to replacing Kij by Pij which evaluates to 1 on the spin 
symmetric state \Q,). Up to terms of the order g^ we have 



7v{dh\n) 

7v{didi+i)\n) 

7T{didi+2)\^) 



[-2g' + Oig')]\n), 
[g^-2g' + 0{g^)]\n) 
[2/ + 0(/)]|0), 



(23) 
(24) 
(25) 



^The same property was postulated [IS] for theta derivative proposed by N. Gromov and P. Vieira 
[5]; here it is a straightforward consequence of the properties of the DunkI operators. 
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while the other combinations quadratic in the dj's do not contribute at the order 
we are considering. Close to the boundaries i = 1,L the result of the projection is 
affected by the boundary conditions and is in general different from the bulk one. 
In particular, we have 7v(di)\il.) = 0, except at the boundary. Collecting the (bulk) 
g^ terms we obtain 



Vl 



1 + 5^^2,2 + 0(/) , 1)2,2 = ^ E (^^ - ^'-+1) ' ' (26) 



with di = dg-. Since we are not particularly interested here in specifying the 
boundary terms, we do not specify the limits of summation above. Let us however 
mention that these terms are important for insuring the properties of the theta 
operator defined in (fT6j) . 

We denote by Dm,n the term appearing at order g"^ and containing nth powers 
of the derivatives. The expression ()26p in the bulk coincides with the result of 
Gromov and Vieira [8] for the B{u) element of the monodromy matrix at two 
loops. At order g'^ we have some contributions from the terms cubic in the Dunkl 
operators 



d^di+i = 


-2/ + 0(/) 


dU,+2 = g' + 0(/) 


(27) 


didf+i -- 


= 2/ + 0(/) 


d^dU = -/ + 0(/) , 


(28) 



where we should understand these relation after projection and action on the vac- 
uum state (that we omit in the following to avoid cumbersome formulae). The 
quadratic and cubic part in the derivatives at order g'^ are reproduced by the 
operator 



A,2 + A,3=E 



3 "^tp, a_ \3 



{d, - d,+2Y - id, - d.+iY + -{d, - d,+2r - 3(9, - dj+l) 



(29) 

This contribution seems rather innocuous, since it vanishes on functions which are 
symmetric in 9. It is absent from j8], where it is postulated that the odd-order 
derivatives are absent. Ruling out the odd-order derivatives might still preserve the 
relation (jlSp . but it will probably introduce unnecessary cross-terms. According 
to N. Gromov, D42 rnay be retrieved in their formulation if the Hamiltonian is 
modified by a transformation which does not change the spectrum. 

For the terms quartic in the Dunkl operators, we have two kind of contributions: 
one is given by two clusters of operators of the type df or didi+i situated in generic 
position with i < j 

djd] = 4/ + 0(/) , d^.ididjdj+i = / + 0(/) , (30) 

djd,dj+i = -2g^ + 0{g^) , did,+id] = -2g^ + 0(/) , (31) 

and the second type is when the two clusters come close to each other (here we list 
only the cases which are different from the general case above) 



dt = Qg^ + 0{g^) 


d-fdU = 3g^ + 0{g') 


(32) 


d^d,+,di+2 = -9' + 0{g^) 


dld,+i = -2,g^ + 0{g^) 


(33) 


dd}+A+2 = 0{g') 


d,df^, = -3g'' + 0{g^). 


(34) 



These contributions can be summarized in 

Z?4,4 = I Y.(d, - d^+l)\d, - a, + i)2 - i Y^id, - a,+i)2(5,+i - d,+2f (35) 



4 
SO that, up to terms of order g^ the operator Vg in the bulk is given by 

Pg^'"^""^ = 1 + g^D2,2 + /(A,2 + 1^4,3 + 1^4,4) + 0(/) . (36) 

To compute the action of Vg on the basis of symmetric sums Pm ■ ■ ■ Pn^ with 
Pn = J2j(^7 it is enough to use the property (fT8|) 

k 

T^e Pni--- PnJe=o = Jl i'^e P^j I^=o) • '^^'^^ 

By direct computation we get Vgpo = L, 'DgP2 = 2g'^L, Vgp^ = Gg'^L. When 
checking the above factorization exphcitly we get some subleading 1/L corrections 
in the products, but they are an artifact of working with finite chain and imposing 
periodic boundary conditions cIl+i = di, which are obviously not compatible with 
the definitions of the Dunkl operators ([!]). As a consequence of (j37p , on an arbitrary 
power of the resolvent Ge{u) defined as 



j 
we have 



] -^ n>0 



Vi GUu] 



d 



e.o V^^^^^(^'^)) • ^''^ 



This relation implies that the substitution G0{u) — )• L-^lnf{u,g) can be done on 
any function which depends on 6 only through the resolvent G0{u). The eigenvalues 
a{u) and d{u) of the diagonal elements A{u) and B{u) on the vacuum can be 
similarly evaluated 

aiu) = /(n+)^ , diu) = f{u-f , (40) 

with u = uibi/2. The function f{u,g), determined above has the same expansion 
to order (7^ as the Zhukovsky variable 

x(^) = i(u+\AI^^4^)=/(^,5) + 0(/)- (41) 

At the order g^, the dispersion relation for the Inozemtsev model starts to be 
different of that of the J\f = 4 dilatation operator, and from that of the BDS model. 
For the full Inozemtsev model, the all-loop expression of the function f{u,g{n)) is 
given [20j by the parametric equation 

/f"")-^'. ^i,>^JE-cJ^]-±iJ^]. (42) 



f{u ) 2TTiK \2kJ 2iK \2k 



with the relation between k and g from eq. ([2]). Here Cii^) is the elhptic zeta 
function with periods 1 and itt/k. 

We beheve that the construction above works for a large class of functions f{u). 
In [25] a procedure was given to construct very general spin chains with long range 
interaction, and the Inozemtsev model falls in the category of chains constructed 
from boost charges in that work. For these models at least, the Dunkl operators 
and the theta derivatives should exist, and we should have 



p„ \n) = {-iTLCn \^) 



with Cn the coefficients of the expansion 



s'"^c' 



1 V^ ^n 

u ^—^ n" 

n>0 



(43) 



(44) 



We hope it is possible to prove this formula by giving an explicit construction of 
the Dunkl operators ([7]). 

The BDS spin chain: The BDS spin chain, defined through its Bethe ansatz 



XiU, 



xiu^) 



M 

n 



Uk - Uj +1 

Uk — Uj — i 



(45) 



was devised such as to reproduce the all-loop dispersion relation of the dilatation 
operator in the su{2) sector. When supplemented with the BES dressing phase 
[35j it reproduces the su{2) dilatation operator at all loop. It was noticed in [21] 
that the BDS ansatz can be formally derived from the inhomogeneous XXX ansatz 
with inhomogeneities given by 9j = 2g sin2ir j/L. In |24j it was shown that the 
BDS spin chain is in fact a projection of the one-dimensional Hubbard model at 
half filling to the spin sector. In that construction, 6j = 2g sin Qj ~ idj are in fact 
dynamical variables corresponding to the momenta Qj of the underlying fermions, 
and in principle subject to backreaction from the spin degrees of freedom. The 
action of the theta operator at all-loop is given by 



n 



BDS 



Pr 



e=o 



i^Pr. 



LC, 



BDS 



with 



C^^' = i2grf^^sin-g, =, Cf^ 



BDs_(25)'"r(n + l/2) 



^^F^(n + 1) 






BDS 

2n+l 



(46) 



0. (47) 



The backreaction decouples from the computation of the symmetric sums, as it 
can be seen in appendix E of [24J. Now, in virtue of the equation ()39p we can 
substitute for any F[G0{u)] 



pr' nceiu)] 



F 



e=o 



L 



^/u^-4.g^_ 
which is of course consistent with the BDS equations ()45p . 



(48) 
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Arbitrary number of magnons: The relations above were given for one- 
magnon eigenstates, but generalization to arbitrary magnon eigenstates is straight- 
forward, due to the morphism property ([20]) . Let us consider the product of M 
copies of the monodromy matrix with auxiliary spaces ai, . . . ,aM- Since the un- 
projected monodromy matrix preserves the space of functions which are symmetric 
under simultaneous permutations of coordinates and spins [27], the product of pro- 
jections is equal to the projection of the product. Proceeding exactly as before, we 
get that 

Ta^{ui) . . .TaJuM)\^) =Ve Ta^,a{ui-9) . . .Ta,,fi{uM]e) \n) . (49) 

From this expression, one can extract any monomial in the operators A{u), B{u), 
C{u), D{u) by appropriate projections in the auxiliary spaces a^. In particular, we 
have 

B{ui)...B{uM)m=VeB^{ui-e)...BQ{uM;9)\^) = \{u})g , (50) 

{n\c{vi)...c{vM) = {^\mc{vi-e)...c{vM;e) = g{{v]\. (51) 

Due to the fact that the R matrix in the Yang-Baxter equation is the same as 
in the XXX model, the algebra of the A{u)^ B{u), C{u), D[u) operators is the 
samcl, with the only modification that the eigenvalues of A{u) and D{u) on the 
vacuum \Q) are now a{u) = f{u^)^ and b{u) = f{u~)^. We conclude that the 
vectors defined in ()50p are eigenvectors of the BDS Hamiltonian if the rapidities 
{u} satisfy the BDS ansatz equations (|45p . 

Scalar products and the correlation functions: In [5] the setup was given 
to compute the correlation functions of the single-trace operators of the AA = 4 
gauge theory in the weak coupling limit. The building blocks entering the correla- 
tion functions are the Slavnov-type scalar products [15j 

S{u;9Uv;e} = {{v;O}\W,0}) (52) 

where \{u;9}) = BQ{ui;9)...BQ{uM',d)\^) and, for example, the rapidities {u} 
obey Bethe ansatz equations with impurities 9. The variables {v} may obey a 
different Bethe ansatz. The tree-level correlation functions can be obtained [5] 
by setting 6 = 0. The scalar products are given by the Slavnov determinant 
formula [151 133 [E] . As observed by direct computation up to two loop order in [H] 
and shown above, at higher loop order the Bethe ansatz eigenstates are deformed 
according to ([50]) into 

IM) ^ \{u})g=Ve\{u;e}) (53) 

and their scalar products can be evaluated using (|39]) . or equivalently (J50l5ip . as 

<?(WIW>9 = Ve{n\Co{vv,e)...Co{vM;0)Bo{uv,e)...Bo{uM;e)\n) 



Vf)S 



e '^{u;e},{v;0} 



e=o 



(54) 



See for example 
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Let us emphasize that, due to the morphism property (j20|l of 'Sq and to the property 
()22p . there are no cross terms in the scalar product, as there are in [S]. This is 
an important difference between our definitions and theirs, and it may give a 
considerable computational advantage at higher loop. The Slavnov determinant 
depends on the impurities d only via the resolvents Gg{u), Go{v). If the theta 
operator is given by the all- loop expression (j46p . the only effect of acting with it 
on the Slavnov determinant is to substitute 

Ge{u) — > L — (55) 

au 

and to insure that {u} (and possibly {v}) obey the BDS ansatz (|45|) . This would 
work even on expressions with finite number of magnons, up to 1/L termq^l- In 
fact, to compute the scalar products we do not really need to know about the 
theta operator; it is enough to know that the monodromy matrix exists and that 
its matrix elements obey the same algebra as for XXX, with modified functions 
a{u) and d{u). This is insured by the morphism properties of the BGHP projec- 
tion. Having an explicit value for the Dunkl operators allows to determine the 
microscopic realization of the model (i.e. to determine the Hamiltonians). 

The quasiclassical limit of the Slavnov determinants, when the distribution of 
mag-nons {u} and {v} condense on some contours Fu and Fv was computed very 
recently by I. Kostov |T7], by using a fermionic representation. The result, which 
generalizes previous results E] is 

ln5{„},{,} = /^Li2(e^(^"(^)+^-(^)-^''(^))-/^Li2(e^^-(^)-^^"(^)) 

Gv Gu 



with Gu,v{z) = / dz' 



(56) 



and the closed contours Cu,v encircling Fu v counterclockwise. We refer to [T7] for 
a discussion of the subtleties in closing the contour Cu around the singularities of 
the dilogarithm. Here we have omitted for simplicity some terms coming from the 
normalization of the states and which can be reconstituted by comparison with 
|17j . These terms will drop out anyway from the physically meaningful quantities 
which are the scalar products divided by the norms. The results at tree level is 
obtained by setting all the impurities to zero. It is straightforward to generalize 
this result at all loop for the BDS model, as well as for the other similar long-range 
models, via the substitution (j55|) . In these cases, we define the quasi momentum 
p(m) through 



g2ipu(M) 



fin 



L M 

^^ U — Uj — I 



In the classical limit u ^ L, so Puiu) is related to the resolvent via 



,„(,) . 0„(,) _ |5^ . (58) 



^ A careful treatment of the boundary terms for the periodic system and of the Hamiltonian insertion 
at the splitting point can be found in [5]. In this reference it was found at two loops, or one loop for 
the correlator, that the two corrections compensate each other. 
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and the Bethe ansatz equations (|45p become in this hmit 

,/ s ^ , -. L din fiz) „ „ , , 

i){z) = ^n , or 0'u(z) = vrn + ^^ , for z£V^. (59) 

Therefore, the results of [17J can be expressed uniquely in terms of the quasi 
momentum and the potential din f (z) / dz . The norm eigenvectors of the Hubbard 
model was conjectured by Gohman and Korepin [Mj; it would be interesting to 
compare the the result at half filling with the prediction for the BDS case. 

Consequence for the correlation functions for AdS/CFT. The all-loop 
Bethe ansatz for the dilatation operator in the su{2) sector can be obtained from 
([15|) by supplementing it with the BES dressing phase (Tbes |3S] • In terms of quasi 
momentum, this amounts to substituting 

p{u) -^ p^^S (^) ^ ^BDS (^) _ • 1j^ ^BES (^) _ (go) 

and in the continuum limit 

pBES(^) ^ ^BES(^) _ ^ . (61) 

2 Vz2 _ 4^2 

When g is large the quasi momentum p^^^{z) reproduces the algebraic curve data 
|36j . It is remarkable that ()56p contains only the quasi momentum p{z) and the 
potential x'{z) = Xj \J z^ — Ag'^ via G0{z). The simplicity of this result leads us, see 
also [17) . to conjecture that the dressing phase can be incorporated in the results 
of |17j simply by substituting the resolvents 

G^,^{z)^G^%^{z). (62) 

According to this conjecture, and assuming that the cancelation between the con- 
tact terms and the Hamiltonian insertions observed in [S] survives at higher loop, 
the all-loop three-point correlators in the su{2) sector would be given by 

lnCu3ig) = -\ E /^Li2[e2^P^(^)] (63) 



1=M, 11,10 ^ 



J{L2-Ll)/2^/z2-4g2+^p^{z) 



J 27r J 2tt 

cs°uCv <^™ 

with the quasimometa Pj{z) defined as 

p^{z) = Gf^''{z)--^^L==, and (u, v, w) ^ (1, 2,3) . (64) 

2 Vz2 - 4s-2 

A general strategy to compute the higher loop correlators was also suggested in 
the conclusions of [5j , where the authors noticed that the substitution of the full 
quasi momentum does not work for finite size systems. Here we conjecture that 
()63p is true in the quasi classical limit, when at least the length Li is large and 
the corresponding rapidities are large, Uj ~ Li, and condense on some cuts. If the 
quasi classical limit of the full solution of the three-point function depends only 
on the quasi momenta, then it should be relatively easy to generalize the result to 
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correlators of three operators in arbitrary positions by using the algebraic curve 
data [Sni ETJ • The recent results [TT] about the Wilson loops and correlators of a 
Wilson loop and a local operator also seem to point out in this direction. The way 
the quasi momentum appears in the classical limit of the Slavnov products reminds 
the thermodynamical Bethe ansatz for the spectrum [38]. It would be instructive 
to clarify the interrelation between the correlators, the amplitudes/Wilson lines 
and the TBA for the spectrum. 

Acknowledgements: The author thanks N. Gromov, O. Foda, I. Kostov, V. 
Pasquier, P. Vieira and D. Volin for illuminating discussions, D. Volin for writ- 
ing a Mathematica code for the Dunkl operators algebra, and N. Gromov and P. 
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